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Abstract 



The neutron-deuteron (nd) doublet ( 2 a n d) and quartet ( a n d) scattering 
lengths have been calculated based on the nucleon-nucleon (NN) interactions 
CD Bonn 2000, AV18, Nijm I, II and 93 alone and in selected combinations 
with the Tucson-Melbourne (TM), a modified version thereof, TM99, and the 
Urbana IX three-nucleon (3N) forces. For each NN and 3N force combination 
also the S H binding energy was calculated. In case of TM99 and Urbana 
IX the 3NF parameters were adjusted to the 3 H binding energy. In no case 
(using np-nn forces) the experimental value of 2 a n d was reached. We also 
studied the effect of the electromagnetic interactions in the form introduced 
in AV18. Switching them off for the various nuclear force models leads to 
shifts of up to +0.04 fm for 2 a n d, which is significant for present day stan- 
dards. The electromagnetic effects have also a noticeable effect on 4 a nc i, which 
otherwise is extremely stable under the exchange of the nuclear forces. Only 
when the electromagnetic interactions are included the current nuclear forces 
describe the experimental value. As a consequence of the failure to reproduce 
2 a n d also the newly measured coherent nd scattering length (b n( i) can not be 
reproduced. For the current nuclear force models there is a strong scatter 
of the 3 H binding energy and the 2 a n ^ values around an averaged straight 
line (Phillips line). This allows to use 2 a n d and the 3 H binding energy as 
independent low energy observables. 
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I. INTRODUCTION 



It has been observed a long time ago that the nd scattering length for total three-nucleon 
spin S=l/2 ( 2 a nd ) is correlated to the 3 if binding energy (Es H ). This correlation is known 
as the Phillips line [1] . Indeed, calculations years later based on simplistic or more realistic 
NN model forces (see [2-7]) yielded quite a few results for the 3 H binding energy and the 
2 a n d scattering length, which lie on or rather close to a line in the two-dimensional plane 
spanned by Ez H and 2 a nd . Also 3N forces of the 27r-exchange type have been added. In [3] 
it was found that this line passes well through the experimental point. 

In recent years chiral perturbation theory and effective theories have been applied to 
nuclear physics. In the pionless formulation [8-10], adequate for extreme low energy phe- 
nomena, it has been shown that 3 H can be energetically stabilized only if a 3N contact 
force is introduced (see however [11,12]). In the two lowest orders of that framework there 
is just one parameter connected to that 3N force. Thus both quantities, Ez H and 2 a n d, 
depend on that one parameter and are therefore correlated though the line does not hit the 
experimental point. In higher orders additional parameters show up and the correlation is 
broken, which makes the two quantities independent. The same observation was made in 
an approach based on chiral perturbation theory [13] which includes explicitely the pion 
degrees of freedom. In the next-to-next-to- leading-order (NNLO), 3N forces occur the first 
time and they depend on two parameters. This makes Ez H and 2 a n d independent and the 
Phillips line correlation is broken. In fact in that framework the two experimental values 
are used to fix those two parameters of the 3N force [13]. Thus we find it interesting to 
ask, whether the conventional and high precision NN forces AV18 [14], CD Bonn 2000 [15], 
Nijm I, Nijm II, and Nijm 93 [16] alone or in combination with the two most popular 3N 
force models , Urbana IX [17] and TM99 [18,19] (an updated Tucson-Melbourne 27r-exchange 
3NF [20] modified in view of chiral symmetry) lead to a strict correlation between Ez u and 
2 a n d or whether that Phillips line correlation is also broken. Further we ask whether the NN 
and 3NF combinations adjusted to Ezh (or may be only one of them) also describe 2 a n d- 
One more reason to confront 2 a n d to state-of-the-art calculations is the recent appearance 
of a precision neutron interferometric measurements of the nd coherent scattering length 
{bud) [21]. 

The coherent scattering length b n d depends in addition to 2 a n d also on the second s-wave 
scattering length for the state of total 3N spin S=3/2 , 4 a n d ■ Because of the Pauli principle 
this quantity is supposed not to be sensitive to short range details of nuclear forces. We 
want to investigate also that quantity in the light of modern nuclear forces. 

Then we would like to add two more investigations. Charge-symmetry breaking (CSB) in 
the strong NN forces is mostly pronounced in the states 1 S'o, where the scattering length for 
the neutron-neutron (nn), a nn , and proton-proton (pp), a pp , systems are different. However 
the value for a nn is still under debate [22,23]. Therefore we would like to present results, 
where the nn forces are replaced by the (strong) pp forces. This will provide some insight 
into the magnitudes of the shifts in 2 a n d caused by small changes in a nn . 

The other investigation is due to effects on the scattering lengths and E 3 h caused by 
electromagnetic interactions, mostly due to magnetic moment interactions (MMI). MMI is a 
relativistic effect and including only that specific force is of course inconsistent, since other 
relativistic effects are not taken into account (see for instance [24] and references therein). 
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But it is interesting to see its separate effect on 2 a n d and A a n d (its effect on the binding 
energy of 3 if and 3 He is known and older results have been reconfirmed recently [25]). Here 
we hit some "defects" in current NN force models. The NN potentials CD Bonn 2000, 
Nijm I, II, and Nijm 93 are fitted directly to the NN data without taking electromagnetic 
interactions (EMI) into account (of course the point Coulomb force in case of the pp system 
has been included). Therefore those strong forces include the effects of the MMI's (and 
further electromagnetic corrections). If we want to see then the separate effects of the 
EMI's we have to subtract from those model NN forces the EMI's and compare to results 
without that subtraction. In case of AV18 the strong force plus EMI's have been fitted to 
the data. Thus the force free of EMI's is just the strong AV18 force alone. In this manner by 
comparing results with and without electromagnetic interactions one can see their effects. 

The paper is organised as follows. In Section II the theoretical formulation is briefly 
outlined. The results are given in Section III, and we end with a summary and an outlook 
in Section IV. More technical details are deferred to an appendix. 



II. FORMULATION 

We use the Faddeev scheme. Including a three-nucleon force a convenient basic formu- 
lation for part of the nd — ■> n + n + p break-up amplitude T [26,27] is the integral equation 

T = tPcj) + (1 + tG ) Vi 1] {l + P)<j> + tPG T + (1 + tG )V/ 4 (1) (l + P)G T. (1) 

The driving term contains the NN t-operator t, permutation operators P, the free 3N 
propagator Go and part of the 3N force, v} 1 ^. Any 3N force can be split into three pieces, 
where for instance the first piece is symmetrical under exchange of particles 2 and 3 , the 
second under 3-1 exchange etc. Thus the quantity vjp is the part symmetrical under 2-3 
exchange like the t-operator t, which is supposed to act on the pair 23. Finally, is the initial 
channel state composed of the deuteron state and a momentum eigenstate of the projectile 
neutron. That integral equation can be solved precisely in a partial wave decomposition in 
momentum space. For details see [28,26,29]. 

The operator U for elastic scattering is given in terms of the amplitude T by quadrature 
as follows 

U = PGo V + PT + \/ 4 (1) (l + P)<j) + \/ 4 (1) (l + P)G T. (2) 

We want to solve the integral equation (1) directly at the threshold of nd scattering. 
This is for zero initial relative momentum if of the projectile and will then lead directly 
to the scattering length. For the convenience of the reader we briefly sketch the necessary 
steps [30]. Our partial wave momentum space basis is denoted by \pqa >, where p and q 
are the magnitudes of standard Jacobi momenta and a a string of angular momentum and 
isospin quantum numbers (see [26,28]). Then for "7f in z- direction we define the auxiliary 
amplitude 



U a ,xi{p,q) = — f — ^rn\Im)(j d m d Im\Jm d + m) < pqa\U\<p > (3) 

m,m d J 
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for the projectile nucleoli with orbital angular momentum A (A = 2A + 1) and total angular 
momentum I combined with the deuteron total angular momentum j d = 1 to total 3N 
angular momentum J. Out of that amplitude one obtains the partial wave projected nd 
elastic scattering amplitude 

uIvm = E / p' 2 dp'Mp')u a > d M{p\ <a>). (4) 

V J 

where <pi(p) are the s- and d-wave components of the deuteron and a d contains the deuteron 
quantum numbers. 

Finally, the projectile spin and the deuteron spin can be combined to the total spin S 
and one obtains 

u}^ - E ^(-y-{l ) £ Wh-Y-{1 V}u} F „. (5) 



The S- matrix element is then given as 



Att 



5 a'E',ae = <Ws'E - i—mq {i) x ' A t/ys',A£ ( 6 ) 



leading to the doublet and quartet scattering lengths for q = 

2n i 

2 a n d = m U*i i 

O u 2' u 2 



2tt 

A a nd = — m Ufa 3- (7) 

One also defines a coherent scattering length b nd as 

bnd = mn + md [ (I) 2 a nd + (-) A a nd } (8) 
m d 6 6 

We defer the special form of the Faddeev integral equation (1) at q = to the appendix. 
It is free of singularities and therefore as easily solved as a bound state problem. Also the 
explicit form of the elastic amplitude for q Q = is given there. 



III. RESULTS 

We used the NN forces CD Bonn 2000, AV18, Nijm I, II and Nijm 93 alone or various 
combinations with the three-nucleon forces Urbana IX, the older Tucson-Melbourne force 
and the modified one TM99. When we combine the Urbana IX 3NF with CD Bonn 2000 
the strength of the repulsive part of this 3NF has been reduced by multiplying it with the 
factor 0.812 in order to get the proper E3 H . 

Due to their nonnegligible influence on the nd scattering lengths, we took special care of 
electromagnetic interactions. In the case of the AV18 potential it is clear how to separate the 
strong AV18 force from the electromagnetic parts because both are well defined and added 
together in fitting the total force to the NN data. In case of the np system the EMI's are given 
in Eqs. (11), (12), and (15) of Ref. [14] and in Eq. (16) for the nn system (for the np system 
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we did not include the very small class IV charge asymmetric force oc L • l/2(lf i — 
Also we neglected the energy dependence of the a 1 ). This is different for the CD Bonn 2000 
and Nijmegen potentials, which were fitted directly to the NN data without adding to them 
electromagnetic interactions, with the exception of the point Coulomb force in case of the 
pp system. Therefore to define the strong forces in the particular NN system one needs 
to subtract the corresponding EMI, which we assumed to be given as in ref. [14]. To be 
precise for the np system we subtract from the np CD Bonn and the Nijmegen forces the np 
EMI's as defined above. Similarly for the nn system we subtract from CD Bonn the MMI 
as defined above. Since we also want to see the effect of replacing the strong nn force by 
the strong pp force we have also to define the strong pp CD Bonn and Nijmegen forces. In 
this case we subtract from those potentials the pp EMI's as given in Eqs. (3)-(8) of Ref. [14] 
without the leading 1 in F c (r) from Eq. (10) of Ref. [14], since the point Coulomb force has 
been taken into account for those potentials fitting to the pp data. 

Before we report on our results we give some comments on our numerical accuracy. 
As usual the partial wave decomposition is truncated at a certain total two-body angular 
momentum j max . Fig. 1 documents the convergence of 2 a nd as function of j max . This shows 
that we reached an accuracy of about three digits. This is for CD Bonn [31] alone and valid 
also for the other NN forces alone. Adding a three-nucleon force we were limited to j max = 5 
due to computer resources. Nevertheless, as Fig. 2 documents, the convergence reached 
for 2 a nd is two digits. In case of 4 a nd with NN forces alone we reach 4 digits convergence 
and including a 3N force an accuracy close to that. This is documented in Fig. 3. The 
other numerical ingredients (discretization in the momenta) are safely under control. In 
all calculations we took into account charge dependence of the NN forces using a simple 
"l^pp(nn) + §W rme t° generate t-matrices in isospin t = 1 2N states [32]. The total isospin 
T = 3/2 3N states have been neglected [32]. We checked that their inclusion does not change 
A a nd up to the fifth digit and the change of 2 a nd is of the order of 0.1%. The triton binding 
energies have been obtained using j max = 6. They are accurate to 2 keV. 

As an overview we show all our results for 2 a nd and Ezn in Fig. 4. We see a group of 
results in the right half of the figure based on NN forces alone and another group close to 
the experimental area including 3N forces. We performed several investigations. First we 
take CD Bonn 2000 as it is (fitted to the NN data) and use the np-nn force combinations 
appropriate to the nd system. In this way the EMI's in the np and nn systems are effectively 
included inside the strong forces. In case of AVI 8 we keep all electromagnetic corrections 
as in [14] except the energy dependence of a' (MMI's for the nn system and the MMI's 
together with the one photon Coulomb term Vc\(np) for the np system). The corresponding 
two predictions are shown as stars in Fig. 4. 

Secondly we want to see the effect of replacing the strong nn forces by the strong pp 
forces. The difference between nn and pp strong forces is mostly located in the different 
scattering lengths a nn and a pp (strong) and will therefore give some information how changes 
in a nn will show up in changes of 2 a nd . Since thereby we do not want to change the EMI's 
we keep in case of AV18 the nn MMI. For the CD Bonn and the Nijmegen potentials the 
strong pp potentials are defined as above and the nn MMI (as for AV18) is added. The 
results are shown as the 5 open circles in Fig 4. As seen in case of CD Bonn 2000 and AV18 
the 3 H binding energy is decreased and a nd increased. For the Nijmegen potentials such 
a comparison is not done, since no nn forces have been introduced. The effects on the 3N 
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binding energies are known. These two first investigations provide theoretical predictions for 
the nd scattering lengths and triton binding energy including all electromagnetic interactions 
similarly as is the case for the measured values. 

In the next two investigations we address the effects of the electromagnetic interactions 
themselves by switching them off while generating theoretical predictions. For the AV18 
potential we take just the np-nn and np-pp strong force combinations alone, while in the 
cases of CD Bonn 2000, Nijm I, II, and 93 we only use the corresponding strong forces 
stripped off from the EMI's, as described above. The resulting theoretical predictions are 
shown as pluses and squares in Fig. 4 for the np-nn and np-pp combinations, respectively. 
Again the binding energy is decreased and a nd increased. 

The individual results of these four investigations are summarized in Fig. 4 also as dashed 
(np-nn with EMI's), dotted (np-pp with EMI's), solid (np-nn), and dashed-dotted (np-pp) 
straight lines fitted in a x 2 sense. These lines include also the corresponding results including 
3NF's (see below). We see a small shift of the lines under exchanges of nn versus pp forces, 
but a more significant shift if the electromagnetic forces are switched off. Though the two 
curves (dashed and dotted) for the cases when the electromagnetic forces are added come 
close to the experimental range spanned by the uncertainty in 2 a n d, they miss it clearly. 
When the electromagnetic forces are switched off the np-nn (solid) and np-pp (dashed- 
dotted) lines go through the experimental point well inside the 2 a nd error bar. 

Now we want to regard our results in more detail as displayed in Table I and in the 
inset of Fig. 4. The theory has to be finally compared to the experimental values, which are 
2 a nd = (0.65 ± 0.04) fm [33], A a nd = (6.35 ± 0.02) fm [33], and b nd = (6.669 ± 0.003) fm [21]. 

The results in Table I are grouped into NN force predictions only and selected com- 
binations with the 3N forces TM, TM99 and Urbana IX. For each potential or potential 
combination we show the results for the various scattering lengths and the 3 H binding ener- 
gies. This is given for the np-pp NN forces, with (without) EMI's in the first (second) row. 
For AV18 and CD Bonn 2000 we also show the results for np-nn forces with (third row) 
and without (fourth row) electromagnetic interactions. Note that in case of the np-nn forces 
including EMI's (as described above) the combinations with TM99 and Urbana IX are fitted 
well to the experimental value -8.48 MeV of the 3 H binding energy. This is also the case for 
Nijm I and II, which, however, refers to the np-pp forces. For the older TM 3N force we 
did not perform a precise (re)fit and the results are only included in view of investigating, 
whether a straight line correlation between 2 a nd and E^h exists. A glance at Fig. 4 tells 
that the individual results scatter around the four straight lines. Thus obviously no straight 
line correlation exist (this has been known before, though for some older calculations the 
numerical accuracy might be not sufficient to give a reliable judgement). 

Let us now concentrate on the group of results with 3N forces. These are displayed 
in the inset of Fig. 4. We see four results (stars) for the np-nn forces including TM99 or 
Urbana IX, where the binding energy has been exactly fitted but where the 2 a nd value is 
too small. These are the results achieved under the supposedly most realistic assumptions 
in this paper. If one switches off the electromagnetic interaction (pluses) the binding energy 
increases and interestingly 2 a nd moves to larger values. Regarding all results, the inclusion 
of the electromagnetic force in our studies shows that they cause shifts of up to about 40 
keV less binding energy and of up to about 0.04 fm decrease in 2 a nd . In no case studied 
the experimental value of 2 a nd is reproduced for np-nn and np-pp strong forces combined 



6 



with different 3NF's with exception of np-pp AV18 combined with TM 3NF, for which the 
theoretical prediction lies at the lower limit of the error bar. 

As one learned from the approach in chiral perturbation theory [13], where two pa- 
rameters are needed to fix the short range 3N forces at NNLO and consequently two 3N 
observables to adjust them, one could foresee that the straight lines in Fig. 4 could only 
by accident pass through the experimental region. For the conventional forces used in this 
paper, one can think of additional 3N force diagrams (the most obvious one the n — p ex- 
change) where a sufficient number of parameters would be available to fit both, Es H and 

Going back to Table I we see that 4 a nd sticks always close to the value 6.34 for the np-pp 
and np-nn NN force choices, without or with 3N forces and with EMI's included. This is 
well within the experimental A a n d error bar. Interestingly, the electromagnetic interactions 
increase A a n d in nearly all cases by about 0.02 and the pure strong force predictions lie always 
outside the experimental error bar. 

Finally, one can confront theory to the very precisely known experimental value of the 
coherent scattering length b n d [21]. Clearly the supposedly most realistic dynamics (nn- 
np NN forces plus TM99 or Urbana IX 3NF's) misses that value. As can be seen from 
Table I, when electromagnetic interactions are included the np-pp force combination reaches 
the experimental value in case of the AV18 and CD Bonn 2000 potentials combined with 
Urbana IX and AV18 with TM99. However, this agreement is accidental and caused by the 
corresponding decrease in 3 if binding. 

IV. SUMMARY AND OUTLOOK 

A recently performed precise neutron interferometric measurement of the nd coherent 
neutron scattering length [21] and a planned precision measurement of the doublet nd scat- 
tering length [34] stimulated us to investigate the theoretical predictions of that quantity for 
the high precision NN forces CD Bonn 2000, AV18, Nijm I, II and 93 in combination with 
currently popular 3N force models. These are the modified 27r-exchange Tucson-Melbourne 
(TM99) and the Urbana IX 3N forces. We have chosen several NN and 3NF combinations, 
which are separately adjusted to the 3 if binding energy. For NN forces alone with and 
without EMI's we recovered the approximate correlation between Ez H and 2 a n d, but the 
scatter around a thought straight line (Phillips line) inside the band spanned by the 4 lines 
in Fig. 4 is quite strong. Adding 3N forces shifts the values into the neighbourhood of the 
experimental range of 2 a n d, but misses the experimental value including its error bar in all 
cases, when electromagnetic forces are included. The inset of Fig. 4 clearly shows that for 
equal or nearly equal 3 H binding energies 2 a n d can vary significantly and vice versa. 

Thus one has to conclude that 2 a n d has to be considered as a low energy observable, 
which is independent from the 3 if binding energy. This observation has been found be- 
fore in approaches based on pure effective field theory (pionless formulation) and on chiral 
perturbation theory (including pion degrees of freedom). Thus in future investigations, ad- 
justing both observables, E-a h and 2 a n d, for conventional nuclear forces will require more 
flexibility in the choice of 3N forces. Adding more mechanisms (on top of the 27r-exchange) 
for 3N forces should be no obstacle. This is a step already performed in the effective theory 
approaches [8-10,13]. 
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We also investigated the effects on 2 a n d resulting from electromagnetic interactions given 
in [14]. The effects on 2 a n d and even 4 a n d are noticeable. For 2 a n d including the electro- 
magnetic interactions reduces its value by up to 0.04 fm. It is interesting to note that 4 a n d 
is perfectly stable under all exchanges of nuclear forces studied in this paper but the elec- 
tromagnetic interactions affect its value, though only in the 3'rd digit. However, only when 
EMI's are included the experimental value is reproduced. 

The effects of adding the electromagnetic interactions on the 3 if binding energy are well 
known and can reach shifts of up to 40 keV less binding energy. 

Due to the failure to describe 2 a n d also the recently newly measured coherent scattering 
length b n d cannot be reproduced theoretically. The good reproduction of A a n d by all inter- 
actions and the small error bar of the coherent scattering length suggests that the value of 
the doublet nd scattering length might be somewhat smaller than the presently one, namely 
around 0.63 fm, different from the present experimental value of 2 a n d- This strongly calls 
for a new, more precise measurement. 

Since the scattering lengths are (extreme) low energy observables, it appears that the 
mentioned effective theory approaches are the most adequate ones. Because one works there 
below a certain momentum cut off, which is smaller than the nucleon mass, they allow 
also to incorporate relativistic effects in a well defined and converged manner. Also 3N 
forces appear in those approaches in a well organised manner, according to a certain power 
counting scheme, and are consistent to the NN forces. In other words, one can take into 
account all these subtle effects, relativity, 3N forces, isospin breaking, in a well controlled 
and systematic manner. In conventional approaches on the other hand, which include a lot 
of phenomenological parametrisations and where no momentum cut-off is used, a reliable 
treatment of relativistic effects poses still a problem and the choices of 3N force mechanisms 
are also quite open. Therefore in conventional approaches physically reliable predictions to 
2 a n d will very likely remain a challenge for quite some time. 



V. APPENDIX 

This appendix summarizes various expressions exactly at the nd threshold qo = 0. The 
first part of the driving term in Eq.(l) turns out to be 

3 q 

< pqa\tP\(j) > = 5 Xo ,o < P^i-^l 2 )^ 1 *' > Vaoiq) 



{lm d I m n \JM) (\ 0^m n \I m n )]j (9) 



The quantities with index refer to the initial state. 
The kernel applied on T is given as 



<pqa\tG PT = J q" 2 dq" 



+ 1 E« Q , ^77 2^a" -JJ, < *25 a \1 

dx ^ ^ (10) 

-1 Xq — X 



with 
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x = ~ kl - r - ^ (11) 

and a contains the same quantum numbers as a with the exception of l a replaced by l a i. 

For our notation see [28]. The deuteron binding energy is written as (—k%/m). The 
remaining parts related to can be worked out correspondingly and can be found in [35]. 
Evaluating the elastic scattering amplitude one needs it at q — go ( see Eq.(4)). Therefore 
the point q = q Q = was included. Then Eq.(lO) simplifies to 



< pq = Oa\tG PT = 2m5 Xafi J q" 2 dq" 



(p, q»; 0) Ea" 2 l «>> g^f«"° < \q"q"a"\T 



One ends up with the elastic scattering amplitude at threshold: 

u XT,xi - -—0x,o0y,o9 ad/ad 2^ —J^\p=o—r\p=o 
III lv p p 

x oi"+l 0/'0/"0 f «2j // / ll\ ^ 1 // // il\rr 

+ h',0 2^ 2 9a d ,*" / 9 d( l m9 ) < „9 9 « K 

+ E / /^V,b'){^4 (1) (l + P)4> + ^ 4 (1) (1 + ^)CoT} Qd/iA/ (p', 0). (13) 
I' J 

The geometrical coefficients g^i 2l ' ll>2 arise from the permutation operator P and are 
given by Eq.(A19) in [28]. 
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TABLES 



TABLE I. Doublet and quartet nd scattering lengths 2 a and 4 a together with the coherent 
scattering length b n d for different NN potentials and selected combinations with different 3NF's. 
All calculations have been done with j ma x = 5. The first and second rows give for each potential 
or potential combination the values obtained with np-pp strong potentials with and without EM 
interactions, respectively (see text for explanation). The third and fourth rows for AV18 and 
CD Bonn 2000 potentials and their combinations with 3NF's are the corresponding results when 
the pp strong NN potential is replaced by the nn one (keeping the nn MMI in case that EMI are 
included). The last column shows our 3 H binding energies. We also included in the second column 
the cut-off parameter A for the TM and TM99 forces. 





A / in 


2 a 


4 a 




El LI 






(im) 


(mil 


(im) 


(MeV) 


CD Bonn 2000 




0.976 


6.347 


6.837 


-7.946 






1.011 


fi 324 

\ > . (J — . i 


fi 833 


-7 989 






0.925 


6.347 


6.812 


-8.005 






0.943 


6.324 


6.798 


-8.048 


CD Bonn 2000+TM 


4.795 


0.622 


6.347 


6.661 


-8.419 




4.795 


0.661 


6.324 


6.657 


-8.463 




4.795 


0.570 


6.347 


6.634 


-8.482 




4.795 


0.590 


6.324 


6.622 


-8.528 


CD Bonn 2000+TM99 


4.469 


0.620 


6.347 


6.660 


-8.422 




4.469 


0.658 


6.324 


6.656 


-8.466 




4.469 


0.569 


6.347 


6.634 


-8.482 




4.469 


0.589 


6.324 


6.622 


-8.527 


CD Bonn 2000+Urb 




0.637 


6.347 


6.668 


-8.423 






0.674 


6.324 


6.664 


-8.467 






0.586 


6.347 


6.643 


-8.482 






0.607 


6.325 


6.630 


-8.526 


AV18 




1.304 


6.346 


7.001 


-7.569 






1.319 


6.326 


6.988 


-7.606 






1.248 


6.346 


6.973 


-7.628 






1.263 


6.326 


6.960 


-7.666 


AV18+TM 


5.215 


0.614 


6.346 


6.656 


-8.478 




5.215 


0.633 


6.326 


6.645 


-8.518 




5.215 


0.556 


6.346 


6.627 


-8.545 




5.215 


0.575 


6.326 


6.616 


-8.584 


AV18+TM99 


4.764 


0.645 


6.346 


6.671 


-8.417 




4.764 


0.663 


6.326 


6.660 


-8.457 




4.764 


0.587 


6.346 


6.643 


-8.482 




4.764 


0.606 


6.326 


6.632 


-8.522 


AV18+UrbIX 




0.636 


6.347 


6.667 


-8.418 






0.654 


6.326 


6.656 


-8.458 






0.578 


6.347 


6.638 


-8.484 
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0.597 


6.326 


6.628 


-8.523 


Nijm I 




1.158 


6.342 


6.924 


-7.742 






1.190 


6.321 


6.919 


-7.782 


Nijm I+TM 


5.120 


0.601 


6.342 


6.646 


-8.493 




5.120 


0.638 


6.321 


6.643 


-8.535 


Nijm I+TM99 


4.690 


0.594 


6.342 


6.642 


-8.485 




4.690 


0.629 


6.321 


6.638 


-8.528 


Nijm II 




1.231 


6.345 


6.964 


-7.663 






1.259 


6.325 


6.957 


-7.700 


Nijm II+TM 


5.072 


0.598 


6.345 


6.647 


-8.500 




5.072 


0.630 


6.325 


6.643 


-8.540 


Nijm II+TM99 


4.704 


0.597 


6.345 


6.646 


-8.487 




4 704 


697 


6 395 


6 649 


-8 597 


Nijm 93 




1.196 


6.343 


6.944 


-7.672 






1.225 


6.322 


6.937 


-7.712 


Nijm 93+TM 


5.212 


0.574 


6.343 


6.633 


-8.502 




5.212 


0.608 


6.322 


6.629 


-8.543 
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FIGURES 




FIG. 1. The convergence of the doublet scattering length 2 a n d as a function of the 2N total 
angular momentum j max for the CD Bonn potential. 
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FIG. 2. The same as in Fig. 1 for the CD Bonn potential combined with the TM 3NF. 
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FIG. 3. The convergence of the quartet scattering length 4 a n d as a function of j max for 
CD Bonn potential (solid curve) and its combination with the TM 3NF (dashed curve). 
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FIG. 4. The results for 2 a n d and E3 H from Table I: np-nn forces alone (pluses), np-pp forces 
alone (squares), and np-nn and np-pp forces plus electromagnetic interactions (stars and circles, 
respectively). The four straight lines (Phillips lines) are x 2 -fhs (np-nn: solid, np-pp: dashed-dotted, 
np-nn with EMI's: dashed, np-pp with EMI's: dotted). The lines with EMI's miss the experimental 
error bar for 2 a nc i [33]. The physically interesting domain around the experimental values is shown 
in the inset. 



16 



